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Abstract 

We derive the evolution equations of parton distribution functions ap- 
propriate in different kinematic regions in a unified and simple way using 
the resummation technique. They include the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi equation for large momentum transfer Q, the Balitskii-Fadin- 
Kuraev-Lipatov equation for a small Bjorken variable x, and the Ciafaloni- 
Catani-Fiorani-Marchesini equation which embodies the above two equa- 
tions. The relation among these equations is explored, and possible ap- 
plications of our approach are proposed. 
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I. INTRODUCTION 



Perturbative QCD (PQCD), as a gauge field theory, involves large loga- 
rithms from radiative corrections at each order of the coupling constant a s , 
such as In Q in the kinematic region with a large momentum transfer Q and 
ln(l/x) in the region with a small Bjorken variable x. These logarithms, 
spoiling the perturbative expansion, must be organized. To organize the log- 
arithmic corrections to a parton distribution function, the various evolution 
equations have been derived. For example, the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) equation sums single logarithms In Q for a large 
x, the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation || sums ln(l/x) for 
a small x, and the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) equation 
0], appropriate for both large and small x, unifies the above two equations. 
The conventional derivation of the evolution equations usually requires com- 
plicated diagrammatic analyses. The idea is to locate the region of the loop 
momenta flowing through the rungs (radiative gluons) of a ladder diagram, 
in which leading logarithmic corrections are produced, and to sum the ladder 
diagrams to all orders. For the DGLAP, BFKL, and CCFM equations, the 
important regions are those with the strong transverse momentum ordering, 
the strong rapidity ordering, and the strong angular ordering, respectively. 

In this paper we shall propose an alternative approach to the all-order 
summations of the various large logarithms. This approach is based on the 
Collins-Soper-Sterman resummation technique, which was developed origi- 
nally for the organization of double logarithms In 2 Q M . Recently, we applied 
this technique to some hard QCD processes, such as deep inelastic scatter- 
ing, Drell-Yan production, and inclusive heavy meson decays, and demon- 
strated how to resum the double logarithms contained in parton distribution 
functions into Sudakov form factors ||. It has been shown [§, that the 
resummation technique is equivalent to the Wilson-loop formalism [0, ||. [| 
for the summation of soft logarithms. Here we shall further show that it 
can also deal with the single- logarithm cases, i.e., the evolution equations 
mentioned above. Therefore, the resummation technique indeed has a wide 
application in PQCD. It will be found that our approach is much simpler 
than the conventional one, and provides a unified viewpoint to the evolution 
equations. 

The procedures of resummation are summarized below. The derivative 
of a parton distribution function with respect to Q or x is first related to 
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a new function, which contains a special gluon vertex. This relation, as a 
consequence of the Ward identity, is exact. Factorizing the new function 
into the convolution of the subdiagram involving the special vertex with the 
original parton distribution function, we derive the evolution equation. At 
this stage, the kinematic orderings of radiative gluons are specified, and the 
subdiagram is identified as the corresponding kernels. We show explicitly 
that the lowest-order subdiagram gives the kernel for the leading-logarithm 
summation. 

In the region with both large Q and small x many gluons are radiated in 
scattering processes with small spatial separation among them, and a new 
effect from the annihilation of two gluons into one gluon becomes impor- 
tant. Taking into account this effect, a nonlinear evolution equation, the 
Gribov-Levin-Ryskin (GLR) equation |10| is obtained. Using the resumma- 
tion technique, the annihilation effect is introduced through next-to-leading- 
twist contributions to the subdiagram containing the special vertex, and the 
GLR equation can be derived easily [[U]]. In the present work we shall not 
address this subject, because it is not very relevant. 

We derive the DGLAP, BFKL and CCFM equations in Sect. II, III, 
and IV, respectively, by means of the resummation technique in axial gauge. 
We explain how the subdiagram containing the special vertex reduces to 
the evolution kernels in the different kinematic regions. Section V is the 
conclusion. 



II. THE DGLAP EQUATION 

Consider deep inelastic scattering (DIS) of a hadron with a light-like mo- 
mentum p = p + 5^ + in the large x limit, where x = —q 2 / (2p ■ q) = Q 2 / (2p ■ q) 
is the Bjorken variable with q the momentum transfer to the hadron through 
a virtual photon. According to factorization theorems, collinear divergences 
arising from radiative corrections to DIS are absorbed into a quark distri- 
bution funtion associated with the hadron, £ being the momentum 
fraction of a parton. The argument p + denotes the large logarithms lnp + 
from the collinear divergences, which will be organized. 

In the axial gauge n ■ A = A + = 0, n M = 5^ being a vector on the light 
cone, the gauge invariant distribution function <p is defined, in the modified 
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minimum subtraction scheme, by 



= / ^e-^ + y-^^a\q(y-)^ + q (0)\P^)U + =o , (1) 

where 7 + is a Dirac matrix, and \p, a) denotes the incoming hadron with 
momentum p and spin a. An average over color is understood. To implement 
the resummation technique, we allow n to vary away from the light cone (i.e., 
n 2 ^ 0) temporarily, and the resultant n-dependent distribution function is 
written as 

(n) (^ + ) = / ^V^ (M | g(f) l 7+g(0) | M) | n ^ (2) 

The above formula is similar to the definition for a quark decay function em- 
ployed in [|J. After deriving the evolution equation, we let n u approach <5 M_ , 
and (j)^ coincides with in Eq. ([]]). That is, the arbitrary vector n appears 
only at the intermediate stage of our formalism, and as an auxiliary tool of 
the resummation technique. < - n ^ can be expressed as the convolution of an 
infrared finite function with 0, and thus they contain the same nonpertur- 
bative information. Their relation will be investigated in details elsewhere 
& 

The key step of resummation is to obtain the derivative p + d^ n > / dp + . 
Because of the scale invariance of < - n - ) in the vector n as indicated by the 
gluon propagator, — iA^(/)// 2 , with 

N^ = g^ + n 2 — — , 3 

n ■ I (n-iy 

(j)^ must depend on p + via the ratio (p • n) 2 /n 2 . Hence, we have the chain 
rule relating p + d/dp + to d/dn 0: 

P + ^ in) = -*-v**«> , (4) 
dp + v ■ n dn a 

with v a = 5 a+ a vector along p. The operator d/dn a applies only to the 
gluon propagator, giving || 

'■ N n» = 1 (i»N av + l v N^) . (5) 



dn a n ■ I 
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Since p flows through both quark and gluon lines, while n appears only 
in gluon lines, the analysis is simplified by considering the derivative with 
respect to n, instead of p + . 

The loop momentum I 11 (l v ) carried by the differentiated gluon contracts 
with a vertex in <j)^ n \ which is then replaced by a special vertex 0] 



n 2 v r 



v ■ nn ■ I 



(6) 



This special vertex can be simply read off the combination of Eqs. (f|) and 
(H). The contraction of (l u ) leads to the Ward identities, 

i{/k+ J) , J _ i fz i{Jz+ /) ( ) 

{k + iy { n k 2 k 2 {k + i) 2 ' 1 ' 

for a quark-gluon vertex, and 



-iN a ^(k + l)^ -iN x ^(k) 



(k + l) 2 ^ uX ' k 2 



-iN a ^(k) -iN a ~t(k + l) 
k 2 (k + l) 2 



(8) 



for the triple-gluon vertex T^x- The Ward identity for the four-gluon vertex, 

^ fabe fade ( n 

Xfiva 



iw « f abe f cde (gx»g^-gx«g,u) + f ace f bde (gx,gv«-gx«g»v) 



+f ade f cbe (gxug^-gx^), (9) 

is simpler. The sum of the four contractions from l x , 1%, 1%, and to the four- 
gluon vertex vanishes: The contractions from l x and 1% to the first term of 
^Xuutj cancel each other, because the first term is antisymmetric with respect 
to the interchange of the indices A and \i. Similar cancellation occurs between 
the contractions from and Zj. The contractions from l\ and 1% and from 1% 
and 11 to the second term cancel separately. The contractions from l x and 11 
and from 1% and l\ to the third term also cancel separately. 

Summing all the diagrams with different differentiated gluons, those em- 
bedding the special vertices cancel by pairs, leaving the one in which the 
special vertex moves to the outer end of the quark line |3| . For a fixed x, we 
obtain the formula, 

p + ^ {n \x ) p + ) = 24>{x ) p + ) , (10) 
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shown in Fig. 1(a), where 4>, the new function mentioned in the Introduction, 
contains one special vertex represented by a square. An identical formula to 
Eq. ([TI]) has been derived for a quark decay function in ||. The coefficient 
2 comes from the equality of the new functions with the special vertex on 
either of the two valence quark lines. Equation fllPf) is an exact consequence 
of the Ward identity without any approximation [Q. An approximation will 
be introduced, when relating (ft to (p^ by factorizing out the subdiagram con- 
taining the special vertex, such that Eq. (|H]) becomes a differential equation 

Of 

It is known that factorization holds only in important regions. The im- 
portant regions of the loop momentum flowing through the special vertex 
are soft and hard, since the vector n does not lie on the light cone, and the 
collinear enhancements are suppressed. In the soft and hard regions can 
be factorized into the convolution of the subdiagram containing the special 
vertex with the original distribution function <f£ n \ 



<t>h{x,P A 



ig 2 C F fj, e 



fi'-'v 



N^ v (l) 



(27r) 4 " e v ■ I 
+2m5{l 2 )<f) {n) {x + l + /p + ,p 



-ig 2 C F /j, € 



SG(j) in) {x,p + ) 



d 4 ~ e l N" u (l) 
-v.. 



(2vr) 4 - e " I 2 



SK(f) (n \x,p^ 



x,p 

) 

/ 



v-l 



(11) 

(n) (x,p + ) 
(12) 



where Cp = 4/3 is the color factor, and SK and SG are additive countert- 
erms. The function S , absorbing the soft divergences of the subdiagram, 
corresponds to Fig. 1(b), where the eikonal approximation for the valence 
quark propagator has been made. The eikonal propagator is represented 
by a double line in the figure. The first and second terms in the integral 
are associated with the virtual and real gluon emissions, respectively, where 
<fi( n \x + l + /p + ,p + ) implies that the parton coming out of the hadron car- 
ries the longitudinal momentum xp + + / + in order to radiate a real gluon of 
momentum Z + . The function absorbing the ultraviolet divergences, cor- 
responds to Fig. 1(c), where the subtraction of the second diagram ensures 
that the involved loop momentum is hard. We emphasize that the factoriza- 
tion formulas in Eqs. ([11]) and ( [12] ) are not exact, but hold only up to the 
leading logarithms lnp + . 
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Employing the variable change £ = x + l + /p + , and performing the inte- 
grations in Eqs. (JTT]) and fll2|) straightforwardly, we arrive at 



4>(x,p~* 



4> s (x,p + ) + 4> h (x,p + ) 



d£ K(x,Z,p + ,n)+G(x,Z,p + , f ,) ^ n \^p + ) , (13) 



with 



K 
G 



— 



1 



up 



Cf In 



1 - x/Z) 
typ 



+ In ^—5(1 - x/i) 
H 



8(1- x/t) 



(14) 



where constants of order unity have been dropped, and v = J (v ■ n) 2 /\n 2 \ is 

the gauge factor. Equation confirms our argument that depends on 
p and n through the ratio (p ■ n) 2 /n 2 = (i>p + ) 2 . In the region with x — *■ 1 the 
logarithm ln(£z/p + //x) in G can be replaced by hi{vp + / n). 

We then treat K and G by renormalization group (RG) methods: 



d K 

fi—K 
dii 



-X 



K 



d r 
-li— G . 

djj, 



(15) 



The anomalous dimension of K is defined by Xk = —fidSK/dfi, whose explicit 
expression is not essential here. When solving Eq. (fl5|), we allow the variable 
fi to evolve from the scale of K to the scale of G. The RG solution of K + G 
is given by 



K{x, f , p + , ji) + G(x, f , p + , fx) 



K(x, £, p + , p + ) + G(x, £, p + ,p + ) 
cp + dpi 

1 



/ -^Aif(a s (/i)) 

Jv+ Li 

C F - 



a s (p 



:i - x/i). 



(16) 



where the initial conditions K(x, £,p + ,p + ) and G(x, £,p + ,p + ) do not contain 
large logarithms after choosing \x as p + . Note that the gauge factor v cancels 
between K and G, implying the gauge invariance of the evolution kernel 
K + G. 
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A remark is in order. The source of double logarithms, i.e., the integral 



containing Xk, vanishes as shown in Eq. (16). If the transverse degrees of 
freedom of a parton are included, an extra factor exp(zlr • b), b being the 
conjugate variable of the transverse momentum k F carried by the parton, 
will be associated with the real gluon emission |4], |5|. Then we deal with a 
two-scale problem. Approximating (fi( n \x + l + /p + , k T ,p + ) by 4>^ n \x, kx,p + ), 
the integration over l + can also be performed, and the convolution in Eq. (|13|) 
is simplified to a multiplication in the 6 space, 

0(x, b, p + ) = [K(x, b, /i) + G(x, p + , n)] (x, b, p + ) , (17) 

with 

K(x,b,fj) = ^^C F \n^- . (18) 

7T 6/i 

Hence, K is characterized by the smaller scale 1/6, and fi should be set to 1/6 
in order to minimize the large logarithms. Consequently, 1/6 is substituted 
for the lower bound of /2 in Eq. fli~6|) , and double logarithms exist. This is 
the difference between this work and || . The above discussion indicates that 
the resummation technique is applicable to the single-logarithm as well as 
double-logarithm cases. 

Inserting Eq. ( |16|) into (|T3|), Eq. ([it]) becomes 

P + -^0 (n W + ) = f T P (x/i)<t> {n) ^P + ) , (19) 

ap^ n Jx 4 

with the kernel 

P(2) = C '(T^)7- < 20 > 

Now we make approach 5^ (the light cone), and obtain 

Q^(x, Q) = ^ f §P(x/£)<Kt, Q) > ( 21 ) 

CLQ 7T Jx 4 

where the variable p + has been replaced by Q, and the gauge invariance of 
the distribution function is restored. It is easy to identify P as the splitting 
function P qq in the limit z — > 1, 

PM = Cf ■ (22) 
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Hence, Eq. ( |i~0|) leads to the DGLAP equation (|2~ip, and the diagrams in 
Fig. 1(b) give the splitting function. 

Note that only the term of P qq singular at z — > 1, ie., £ — > x, is reproduced 
in the resummation approach. The reason is as follows. The z 2 term in the 
numerator of P qq arises from the radiative correction with the ends of a real 
gluon attaching each valence quark line, whose loop integral is written as 



2/7 

9 C F 



dH 



7m 



(2vr) 4 

xN^(1)2tt5(1 2 )5 [i-x 



{i P - if 2 (& - if 



p 



(23) 



fp in the trace comes from the Dirac structure of the proton distribution 
function. The numerator £ /p— /I contributes the factor £p + — l + = xp + 
because of the last 5 function. The denominator (£p — I) 2 = —2£p + v ■ I gives 
the factor £p + due to the on-shell gluon with I 2 = 0. Their combination leads 
to a power of z = x/£. The two fermion propagators then explain the term 
z 2 . 



When applying the derivative in Eq. 
first term renders Eq. (|23|) reduce to 



to N^ u , we have Eq. (gj), whose 



-fa 



dH 



-Tr 



(2nf 

xv a N au (l)27r5(l 2 )5 f £ 



2 (ep-Z) 2 
/+ 



X 



P 



(24) 



as described by Fig. 2. To obtain the above expression, we have used 



Zt>- [ 

(£p - 1) 2 



-2&+v ■ I 



-1 + 



2£p+v ■ I 



(25) 



where the second term gives a vanishing contribution when multiplied by jp 
in the trace. Hence, the differentiation with respect to n employed in the 
resummation technique removes a power of z. If further applying the eikonal 
approximation to the remaining fermion propagator, 



lv 



-2£p+v ■ I 



lv 



2jp + v v 
-2&+v ■ I 



v-l 



(26) 
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Fig. 2 reduces to the seccond diagram in Fig. 1(b), and another power of z 
disappears. 

We emphasize that the resummation technique is based on PQCD fac- 
torization theorems, and the factorization of the subdiagram makes sense 
only in the soft and hard regions. That is, the function (f> a , with the eikonal 
approximation, collects the leading contribution from the soft region £ — > x, 
and 4>h collects the leading contribution from the hard region. The subdi- 
agram in Fig. 2, containing a finite piece of the splitting function, can not 
be absorbed into 4> s or <fih, since the quark propagator is not eikonalized, 
and the real gluon is not hard. Therefore, the finite part of the splitting 
function is completely missing in the resummation approach. In the conven- 
tional derivation of the DGLAP equation one-loop diagrams are computed 
explicitly without resort to eikonal (soft) approximation, and this makes the 
difference. The same conclusion applies to other cases that involve a soft 
approximation, such as the conventional derivation of the CCFM equation 
H], where only the terms of the relevant splitting function singular at z — > 1 
and at z — > were reproduced, and the finite part was in fact put in by hand, 
as mentioned in Sect. IV. 



In this section we demonstrate that the resummation technique reduces to 
the BFKL equation for the gluon distribution function in the small x region. 
The unintegrated gluon distribution function F(x, kx), defined by 



in the axial gauge n ■ A = 0, describes the probability of a gluon carrying a 
longitudinal momentum fraction x and transverse momenta k^. is the 
field tensor. Similarly, we vary the vector n arbitrarily first to work out the 
resummation, and then show that the BFKL kernel is independent of n as in 
the DGLAP case. After deriving the evolution equation, n is brought back 
to the light cone. 



III. THE BFKL EQUATION 




x-Y,(p^\^(y',yT)F^(0)\ P ,a) , 



a 



(27) 
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We do not make explicit the p + dependence of F for the following reason. 
The BFKL equation governs the behavior of F with the momentum fraction 
x. The variation of x can be achieved by varying the hadron momentum 
p + , if the parton momentum xp + = k + , appearing only in the exponent 
in Eq. (|27|), is fixed. Again, F involves the vectors k and n, which should 
combine into the ratio (k ■ n) 2 /n 2 . Since k + is fixed, we do not choose it 
as an argument of F, and regard that F depends on p + implicitly through 
x = k + /p + . Hence, the derivative of F with respect to x is related to the 
derivative with respect to p + considered in the resummation formalism, and 
the chain rule in Eq. (0J) holds: 

d d n 2 d 

- x—F(x, k T ) = p + T -F(x, k T ) = v a — F(x, k T ) . (28) 

dx dp + v ■ n dn a 

Applying the operator d/dn to a gluon propagator, we get Eq. (|) and the 
same special vertex in Eq. Summing all the diagrams with different 

differentiated gluons and employing the Ward identity, the special vertex 
moves to the outer end of the parton line as explained in Sect. II j4|. We 
then obtain the derivative of F, 

- x-^F(x, p T ) = 2F{x, p T ) , (29) 

described by Fig. 3(a), where the new function F contains one special vertex. 
It is easy to observe that Eq. (|29|) is the copy of Eq. fllO|) for the gluon 
distribution function. 

We relate F to F by factorizing out the subdiagram containing the special 
vertex, such that Eq. reduces to a differential equation of F. In the 
leading soft and hard regions of I, which flows through the special vertex, 
the factorization is performed according to Figs. 3(b) and 3(c), respectively. 
Fig. 3(b) collects the soft divergences by eikonalizing the gluon propagator. 
We extract the color factor from the relation fabcfadc = —N c 5 a d, where the 
indices a,b, . . . have been indicated in the figure, and N c = 3 is the number 
of colors. The corresponding factorization formula is written as 



F s (x,k T ) = iN c g 2 J— -T^vpl-iN^ 



(2tt) 4 ^ PL V /J -2xp-l 

' (30) 



2 m 5(l 2 )F(x, |k T + 1 T |) + ^ F(x, k T ) 
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where iN c comes from the product of the overall coefficient — % in Eq. (^) 
and the color factor —N c extracted above, and the triple-gluon vertex for 
vanishing I is given by 



(31) 



The first term in the brackets corresponds to the real gluon emission, where 
F{x, |ky + It |) implies that the parton coming out of the hadron carries the 
transverse momenta kr + 1^ in order to radiate a real gluon of momenta 
It- The second term corresponds to the virtual gluon emission, where the 9 
function sets the upper bound of It to kx to ensure a soft momentum flow. 

It can be shown that the contraction of p with a vertex in the quark 
box diagram the partons attach, or with a vertex in the gluon distribution 
function, leads to a contribution down by a power 1/s, s — (p+q) 2 , compared 
to the contribution from the contraction with Vp. Following this observation, 
Eq. flBTf) is expressed as 



F s (x, k T ) 



dH 



iN c g 2 

+ T 2 



{2K)A N^(l) V £f[2ni5{l 2 )F(x,\k T + l T \) 



■F(x, kj 



(32) 



which corresponds to Fig. 3(b) exactly. The eikonal vertex v„ comes from 
the last term divided by xp + in Eq. ([TTJ). The remaining metric tensor g m 
has been absorbed into F. Assuming n = (n + ,n~,0) for convenience, the 
integrations over V~ and / + to infinity in Eq. (|3"2"D give 



F s (x, k T ) 



2 

x 

ah 
2 



d 2 h 



-n 



\l+=oc 
\l+=0 



F(x, |k T + l r |) - 6(k 2 T - l 2 T )F(x, k T ) 



d 2 l T 



F(x,\k T + l T \)-6(kZ-%)F(x,k T ) , (33) 



with a s = N c a s /ir. The first line of the above formulas demonstrates explic- 
itly how the n dependnece cancels in the evaluation of Fig. 3(b). 
The contribution from the first diagram of Fig. 3(c) is written as 



T 



d 2 u 



71 



1 



l\ + (Jfe+i/) 
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k+u 



In 



1% + (k+u) 2 



k+ 



v 



2 [if, + (k+^ ^2 + {k+u) 2 + k+u 



(34) 



which confirms the statement that F depends on x (or p + ) via the ratio 
(k ■ n) 2 /n 2 = (k + u) 2 . In the interesting region with small k + , Eq. ( [54]) is less 
important (it vanishes as fc + — > 0), and _F S dominates. Ignoring the contri- 
bution from Fig. 3(c) along with its soft subtraction (the second diagram), 



that is, adopting F ~ F s , Eq. 



becomes 



dF(x, k T ) 
d\n(l/x) 



a. 



(fir 



F(x, |k T + l r |) - 0(4 - l T )F(x, k T ) 



(35) 



which is exactly the BFKL equation. The n dependence residing in Fig. 3(c) 



is removed with the vanishing of Eq. (34) at small k + , and the BFKL ker- 
nel turns out to be gauge invariant. It is understood that the diagrams in 
Fig. 3(b) play the role of the BFKL kernel, a similar conclusion to that drawn 
at the end of Sect. II. 

Since the explicit dependence of F on the large scale p + is neglected, the 
transverse degrees of freedom of a parton must be taken into account. This 
explains why the gluon distribution function at small x is constructed based 
on the high-energy /^-factorization theorem |13]. Therefore, our formalism 



is applicable to the distribution functions defined according to the collinear 
factorization (the DGLAP case) and according to the /cT-factorization (the 
BFKL case). 



IV. THE CCFM EQUATION 

With the discussion in the previous two sections, it is not difficult to 
demonstrate that the resummation technique reduces to the CCFM equation, 
which embodies the DGLAP equation and the BFKL equation. It hints that 
we should maintain both the l + and It dependences in the unintegrated 
distribution function for the real gluon emission, namely, consider F(x + 
l + /p + , |k T + l T \). The BFKL equation is appropriate for the multi-Regge 
region, where the transverse momenta carried by the rung gluons of a ladder 
diagram are of the same order, i.e., It ~ fcr- Hence, the loop momentum It 
flowing through the distribution function is not negligible, and we can make 



13 



the soft approximation 

F(x + l + /p + , |k T + 1 T |) « F(x, |k T + 1 T |) , (36) 

in F s , from which the BFKL equation is derived. While the DGLAP equation 
is appropriate for the transverse momemtum ordered region, in which we have 
It <C pr, and thus the approximation 

F{x + l+/p+, |k T + 1 T |) « F(x + fey) . (37) 

The kj< dependence of the distribution function decouples, and is integrated 
out from both sides of Eq. (|35l). This is the reason a parton distribution 



function in the DGLAP equation needs not to involve the transverse degrees 
of freedom. The argument x + l + /p + then leads to the splitting function. 
The same diagrams in Figs. 1(b) and 3(b) give the different evolution kernels, 
because they are factorized according to the different kinematic orderings of 
radiative gluons. 
Start with 

p + -^F(x, k T ,p + ) = 2F(x, k T ,p + ) , (38) 

where the arguments hr and p + of the unintegrated gluon distribution func- 
tion manifest the attempt to unify the BFKL and DGLAP equations. Sim- 
ilarly, the new function F involves one special vertex at the outer end of a 
parton line. If following the standard procedures of resummation, we should 
factorize out the subdiagram containing the special vertex in the leading 
soft and hard regions, and derive F s and Fh, respectively, as in the previous 
sections. The function Fh involves the lowest-order virtual gluon emission. 
This idea leads to a new unified evolution equation, which will be studied 
elsewhere [fL4|. To reproduce the CCFM equation, however, the inclusion of 



virtual gluons must be performed in a different way: They are embedded in 
F s , instead of absorbed into Fh- Consequently, the factorization of the subdi- 
agram is described by Fig. 4(a), where the two jet functions J group all-order 
virtual corrections, and the lowest-order real gluon emission between them 
is soft. 

First, we resum the double logarithms contained in J by considering its 
derivative 

P + ^^(Pt,P + ) = J(k T} p + ) 

= [Kj(k T ,fi) + Gj(p + ,fi)]J(k T ,p + ) , (39) 
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which is similar to Eq. (O). At lowest order the function Kj comes from the 
first diagram of Fig. 3(b), and Gj from the two diagrams in Fig. 3(c). The 
relation between Kj + Gj and J is simply multiplicative, since J collects only 
virtual gluons. We have set the infrared cutoff of Kj to k T , as indicated by its 
argument. This cutoff is necessary here due to the lack of the corresponding 
real gluon emission, which serves as a soft regulator. The one-loop Kj can 
be obtained simply by working out the second integral in Eq. (^) without 
the 9 function. The anomalous dimension of Kj is then found to be 7j = a s . 
The function G j can also be computed, but its explicit expression is not 
important. The standard RG analysis leads to 



fP + dfi 

Kj(k T , /i) + Gj(p + , fj) = - / —jj(a s {fl)) 

Jk T /! 



(40) 



with the initial conditions Kj(Ict, hr) = Gj(p + ,p + ) = 0. Of course, we have 
neglected the constants of order unity in Kj and Gj. 
Substituting Eq. ([40]) into (|39|), we solve for 



J(k T ,Q) = A(Q,k T )J {0) 



with the double-logarithm exponential 
A(Q,fc T 



exp 



-a.s 



Q dp + fP + dfi 

k T p + Jk T U 



(41) 



(42) 



We have chosen the upper bound of p + as Q, and ignored the running of a s . 
The initial condition can be regarded as a tree-level gluon propagator, 
and then eikonalized in the evaluation of the soft real gluon emission below. 
We split the above exponential into 



A(Q,fc. 



T 



A 



l/2 (Q,zq)A^ 2 s (z,q,k T ) 



(43) 



with z = x/£ and q = It / '(1 — z) , where £ is the momentum fraction entering 
J from the bottom, and It is the transverse loop momentum carried by the 
real gluon. The so-called "Sudakov" exponential As and the "non-Sudakov" 
exponential A^vs are given by 



As(<2, zq) = exp 



-2a, 



Q dp + fP + d/i 

zq p + Jk T fL 
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A 



NS[Z,q,PT) 



exp 
exp 
exp 



-at. 



-2a, 



-as. 



/•Q 2 dp 2 


rl-hr/p 


dz' 


/(2<Z)2 J) 2 . 




1 


- Z 




fP + dp,' 








Jk T ft 








r k r dp 2 ' 




/2 Z 


J(z'q) 2 P 


2 





(44) 



where the variable changes ft — (1 — and p + = p for As, and ft — p and 
p + = z'g for Atvs have been employed. 

With Eq. Fig. 4(a) reduces to Fig. 4(b), where the tree-level gluon 
propagator on the right-hand side has been eikonalized, and that of the 
left-hand side has been absorbed into F. Based on Fig. 4(b), F is written as 



F(x, kTiP^ 



x6(Q - zq)F{x + l + /p + , |k T + 1 T \,P A 



(45) 



Basically, the above formula is similar to the real gluon emission term in 
Eq. (|32|) except for the exponential A 2 from the two jet functions J, and 
the 9 function. The 9 function requires Q > zq such that the Sudakov 
exponential As is meaningful, which comes from the angular ordering of 
radiative gluons Q/(xp + ) > It/[(£ ~ x)p + ]. Compared with the transverse 
momentum ordering for the DGLAP equation || [15|], Q (It) has been divided 
by the gluon energy xp + ((i — x)p + ). Hence, the inserted scale zq reflects the 
special kinematic ordering for the CCFM equation. Those radiative gluons, 
which do not obey the angular ordering, conntribute to the non-Sudakov 
exponential. This is one of the motivations to introduce the scale zq. 

Using the variable change £ = x + l + /p + and performing the integration 
over l~ , we obtain 



a. 



F(x,k T ,p + ) = — 



di 



d 2 h 



2n 2 (£-x)p +2 



vr [n+l 2 T + 2n-{£-x) 2 p+ 2 ] 2 
x9(Q-zq)F(Z,\k T + \ T \,p + ) , 



A 2 (Q,k T ) 



(46) 



where n = (n + ,n~,0) has been assumed. Equation ( ]46|) is then substituted 
into (^) to find the solution of F. We integrate Eq. (fIBD from p + = to Q, 
and adopt the variable changes £ = x/z and ly = (1 — z)q. To work out the 
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p + integration, F(x/z, |k-r + lT|,P + ) is approximated by F(x/z, \pt + It\,It)- 
The solution of F is given by 



1 A^ 1 

F{x,k T ,Q) = F^ + a s [ dz { -^e{Q-zq)A s {Q,zq)A NS {z,q,k T 

Jx J irq z 



1 



x 



z{\ - z) 



F(x/z,\k T + (l-z)q\,l T ) 



(47) 



where the nonperturbative initial condition F^ 3 ' corresponds to the lower 
bound of p + . Again, the n dependence disappears for a similar reason to 
that for Eq. (P5|), and the CCFM kernel is gauge invariant. 
Equation fllTD can be reexpressed as 



F{x,k T ,Q) 



F (o) + f 1 

Jx 



dz 



<Pq_ 
nq 2 



9(Q-zq)A s (Q,zq)P(z,q, k T ) 



xF(x/z,\k T + {l- z)q\,l T ) 



with 



P 



a. 



1 



l-z 

which is close to the splitting function 



+ A NS (z, q,k T )--2 + z(l - z) 



P, 



a. 



1 1 
+ - 

— z z 



2 + z(l 



(48) 
(49) 

(50) 



Obviously, Eq. @ is the CCFM equation @. To arrive at Eq. (^y), we 
have employed the identity 1/(^(1 — z)) = 1/(1 — z) + 1/z, and put in by 
hand the last term — 2 + z(l — z). This term, finite at z — > and at z — > 1, 
can not be obtained in the conventional approach either || as stated at 
the end of Sect. II. Note that only the non-Sudakov form factor Ajvs in 
front of 1/z is kept, because A NS vanishes when the upper bound zq of p + 
approaches zero, as shown in Eq. (fill) , and thus smears the z — > pole of 
the function P. This is another motivation to introduce the scale zq. The 
above derivation shows that the complicated diagrammatic analysis involved 
in the conventional derivation of the CCFM equation [[J is avoided using the 
resummation technique. 



VI. CONCLUSION 
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In this paper we have shown that the resummation technique provides a 
unified and simple viewpoint to the organization of the various large loga- 
rithms, and reduces to the DGLAP equation, the BFKL equation, and the 
CCFM equation in the different kinematic regions. The main idea is to relate 
the derivative of a parton distribution function to a new function involving 
a special vertex. The summation of the large logarithms is embedded in the 
new function without resort to complicated diagrammatic analyses. When 
expressing the new function as a factorization formula, we obtain the evo- 
lution equation. The subdiagram containing the special vertex, factorized 
according to the specific orderings of radiative gluons, leads to the corre- 
sponding kernels. By means of the resummation technique, the connection 
among the evolution equations becomes transparent. 

Many applications of our formalism to small x physics follow this work, 
which will be briefly described below. In Eq. ( |5D| ) associated with the BFKL 
equation, the real gluon emission term in fact involves the distribution func- 
tion F(x + l + /p + , |k T + 1 T |) as in the DGLAP and CCFM cases, instead of 
F(x, |ky + lx|). The former is replaced by the latter when the strong rapidity 
ordering x + l + /p + ^> x i.e., Eq. ( j36|) is employed. Taking this into account, 
it is of no doubt that the loop momentum / + should not be extended to 
infinity in the derivation of Eq. (0), since F{x + l + /p + ) vanishes at large 
momentum fraction. If truncating l + at the scale of order Q, the resultant 
BFKL equation contains an intrinsic Q dependence. This modified equation 
has been proposed and studied in [|16| . The HERA data of the DIS structure 
function F 2 (x,Q 2 ) which exhibit a steep rise (corresponding to hard 



pomeron exchanges) at small x for large Q 2 ~ 50 GeV 2 , and a flat rise for 
low Q 2 ~ 4 GeV 2 (because soft pomeron contributionns begin to play), were 
explained successfully. Note that the Q dependence of the data can not be 
understood by means of the conventional BFKL equation. 

Another consequence of relaxing the rapidity ordering is the recovery of 
the unitarity of the BFKL evolution. It can be easily recognized that real 
gluon emissions are responsible for the rise of the gluon distribution and the 
structure function F 2 - However, the approximation in Eq. fl5E|) overestimates 
the real gluon contributions, such that F 2 rises as a power x~ x , and violates 
the unitarity bound F 2 < const. In 2 (l/x). By adopting F(x + l + /p + ) in the 
evaluation of the BFKL kernel, we have been able to show that the power rise 



is moderated into a logarithmic rise The reason we have more freedom 



to modify the BFKL equation is that the gluons are not reggeized according 
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to the rapidity ordering before the BFKL kernel is derived. 

The BFKL equation for the polarized gluon distribution function can 
be derived simply in our formalism |Tj|. After factorizing the subdiagram 



containing the special vertex, we contract Eq. (|30|) with the gluon polarization 
vectors e M e 7 for positive and negative helicities, and take their difference. 
The obtained equation will be employed to study the small x behavior of the 
polarized structure function gi(x), and may help to clarify the discrepancies 
between the previous and recent data of g% at x ~ 1CT 2 , and between the 



Regge extrapolation and the QCD fit for gi at x — > pO . 

It is well-known that the soft divergences from virtual and real radiative 
corrections to DIS cancel each other order by order. However, in the deriva- 
tion of the CCFM equation virtual gluons are summed to all orders and 
grouped into the Sudakov and non-Sudakov form factors, while real gluon 
contributions are evaluated only to lowest order, which leads to the split- 
ting function P gg . To fulfill the soft cancellation, the coupling constant a s 
in the exponent of the Sudakov form factor As must be frozen, such that 
the lower bound q of the variable p + can reach zero, i.e., As — > as shown 
in Eq. (fHj). The infrared singularity from the factor 1/q 2 in the kernel of 
Eq. ( f48|) is thus suppressed. Otherwise, a running a s (fl) will prohibit q from 
being below Aqcd- However, the constant a s becomes a parameter. With 
our approach, it is trivial to derive a new unified evolution equation for the 
InQ and ln(l/x) summations, in which both virtual and real corrections are 



considered to lowest order 14 



At last, to improve the accuracy of the kernel to next-to- leading loga- 
rithms, we only need to evaluate the 0{a 2 s ) subdiagram. Such an evaluation 
can be performed in a straightforward way. The BFKL equation including 
the summation of the next-to-leading ln(l/x) will be published elsewhere. 

This work is supported by National Science Council of R.O.C. under the 
Grant No. NSC87-2112-M-006-018. 
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Figure Captions 



FIG. 1. (a) The derivative p + d(f)/dp + in the axial gauge, (b) The soft struc- 
ture and (c) the ultraviolet structure of the 0(a s ) subdiagram containing the 
special vertex. 

FIG. 2. One of the diagrams that contributes to the finite part of P qq . 

FIG. 3. (a) The derivative —xdF/dx in the axial gauge, (b) The soft struc- 
ture and (c) the ultraviolet structure of the 0(a s ) subdiagram containing the 
special vertex. 

FIG. 4. (a) The subdiagram containing the special vertex for the CCFM 
equation, (b) The subdiagram for the CCFM equation after resumming the 
double logarithms in J. 
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